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Abstract 

The general theory of motion in the vicinity of a moving quantum nodal point (vortex) is studied 
in the framework of the de Broglie - Bohm trajectory method of quantum mechanics. Using 
an adiabatic approximation, we find that near any nodal point of an arbitrary wavefunction ijj 
there is an unstable point (called X-point) in a frame of reference moving with the nodal point. 
The local phase portrait forms always a characteristic pattern called the 'nodal point - X-point 



complex'. We find general formulae for this complex as well as necessary and sufficient conditions of 



validity of the adiabatic approximation. We demonstrate that chaos emerges from the consecutive 
scattering events of the orbits with nodal point - X-point complexes. The scattering events are of 
two types (called type I and type II). A theoretical model is constructed yielding the local value 
of the Lyapunov characteristic number in scattering events of both types. The local Lyapunov 
. characteristic number scales as an inverse power of the speed of the nodal point in the rest frame, 



implying that it scales proportionally to the size of the nodal point X- point complex. It is also 



' an inverse power of the distance of a trajectory from the X-point's stable manifold far from the 

' complex. This distance plays the role of an effective 'impact parameter'. The results of detailed 

numerical experiments with different wavefunctions, possessing one, two, or three moving nodal 
^ ' points, are reported. Examples are given of regular and chaotic trajectories, and the statistics 

of the Lyapunov characteristic number of the orbits are found and compared to the number of 
encounter events of each orbit with the nodal point X-point complexes. The numerical results are 
in agreement with the theory, and various phenomena appearing at first as counter-intuitive find 
a straightforward explanation. 
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I. INTRODUCTION 



The dynamics in quantum systems with vortices, i.e. singularities of the phase field 
of the wavefunction tjj = i?e ' , has attracted much interest in recent years be- 



cause oi 
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a variety of potential applications, e.g. in tunneling through potential barriers, 
aUistic electron transport 5,0,1^1, superfluidity 10 1, Bose- Einstein condensates 
14. Il5|. optical lattices [l6|, atom-surface scattering ITJ, Josephson junctions 
li, decohe.enc:H et. The so-called •t.aiecto.y approach lend, itself very co,„e„,e„tly 
to such a study (see [20| for a review). In this approach one follows the orbits of 'par- 
ticles' tracing the quantum-mechanical currents. This is computationally equivalent to a 
Lagran gian quantum-hydrodynamical approach 2l| or to the 'Bohmian' or 'pilot wave' ap- 
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2J]. The trajectories are described by first order equations of the form 



proach 

V = VS{ x)/h. Such trajectories provide a Lagrangian visualization of quantum processes 



[e.g. 



25| ) which is distinct from the Eulerian (i.e. Schrodinger) approach, although it is 



consistent with it. 

A number of authors have found that the quantum trajectories of low-dimension al sy stems 



can be very chaotic [26 
physical importance o 
recent papers of ours 
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4^. The 



chaotic quantum trajectories has been extensively discussed in three 
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The generation of chaos is directly associated with the appearance of quantum vortices. 
It has been pointed out 39|, |4lJ that chaos is, in general, caused by the motion of quantum 
vortices. In the case of fixed vortices, on the other hand, chaos can still be generated if we 



allow a coupling of the wavefunction to a vector electromagnetic potential (e.g. 3J, l46l|). 



In previous papers of ours ([44|, hereafter EKC, and [45|) we studied the quantum trajec- 
tories in particular examples of 2D systems with a moving quantum vortex, in order to find 
numerical indications of the mechanism of generation of chaos. We note that this problem is 
quite different from traditional problems of nonlinear dynamics. First, the equations of mo- 
tion become singular near a vortex. Furthermore, the vortex is oscillating quasi-periodically, 
i.e., with more than one incommensurate frequency. Thus there are no obviously identifiable 
critical points of the flow other than the vortex itself. 

Setting the identification of the critical points as a primary target, in EKC we looked for 
such points in a moving local frame of reference centered at a moving nodal point, and made 
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use of an approximative form of the equations of motion valid under a so-called adiabatic 
approximation. Our main finding can be summarized as follows: In the above frame and 
approximation, the nodal point is seen to create a saddle point nearby, called the 'X-point'. 
The local phase portrait was called the 'nodal point - X-point complex'. Most trajectories do 
not penetrate deeply inside the complex. Instead, they are accelerated along the X-point's 
asymptotic curves and eventually they are scattered by the complex. Chaos is generated 
by a sequence of such scattering events. This conclusion was substantiated by following 
the evolution of the deviation vectors of some representative chaotic trajectories. We also 
found the domains avoided by the nodal point - X point complex and demonstrated that the 
trajectories covering these domains are regular, i.e. they obey effective integrals of motion 
and they have zero Lyapunov characteristic numbers. 

Our study so far relied only on numerical examples of trajectories in particular -0— fields, 
in which (i) the wavefunction ip{x,y,t) has a simple form, and (ii) there is only one nodal 
point present in the configuration space {x, y) at any time t. These restrictions are removed 
in the present paper, in which (i) we develop the general theory of motion near moving 
2D quantum vortices, applicable to generic ip— fields, and (ii) we continue the study of 
particular numerical examples, in cases with more than one nodal point. 

Regarding (i), section II contains a general analytical treatment of the motion near the 
critical points of the quantum flow yielding: a) the form of the equations of motion in terms 
of the coefficients of a local expansion of a generic -0— field around a nodal point, b) general 
formulae for the structure and stability of the nodal point - X-point complex, c) conditions 
of validity of the adiabatic approximation, and (most importantly) d) theoretical predictions 
for the values of the Lyapunov characteristic numbers generated locally by the interaction 
of the trajectories with a nodal point - X-point complex. The latter problem is treated like 
a classical scattering problem. The most important parameters in the theory turn to be the 
speed of the nodal point and the impact parameter, i.e. distance of a trajectory from the 
X-point's stable manifold far from the complex. The theory leads to a quantification of the 
degree of chaos, i.e. the level of the Lyapunov characteristic number of the trajectories, in 
systems with quantum vortices. 

Regarding (ii). Section III tests the theory of section II against numerical experiments, 
focusing on examples in which the -0— field generates more than one nodal point at the same 
time. We obtain numerical values of the Lyapunov characteristic numbers for statistical 
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ensembles of orbits and compare these values with the predictions of section II. We also check 
the quantitative relation between the Lyapunov characteristic number and the number of 
encounters of a trajectory with the nodal point - X-point complexes. Section IV summarizes 
our conclusions. 



II. THE MOTION IN THE VICINITY OF A NODAL POINT 



A. Equations of motion 



Let ^a;o(t), Voit) j represent the center of a moving frame of reference, V{t) = Vy) = 
(xo, i/o) being its velocity at the time t with respect to the rest frame (x, y). The wavefunction 
can be expanded around (xo,?/o)- Up to second degree we have: 

+ ^(^ao2{t)+ibo2it)^v^ + (^anit)+ibnit)^uv + ... (1) 

where u = x — Xq, v = y — i/q and the coefficients aij, bij are real. 

The equations of motion in such a frame, which follow from the equations of motion 
X = VS {h = 1) in the rest frame, read: 

(^»=/m(^)-(K.,K.) . (2) 

All the frames of reference moving with the same velocities Vx{t), Vy{t) at all times t form 
an equivalence class with respect to parallel translations in the configuration space. We 
consider as representative of the class a frame, of which the center ^a;o(to)) 2/o(^o)^ coincides 
at some time t = to with the instantaneous position of a nodal point {xj\f(tQ),yj^(tQ)) of the 
wavefunction ip, i.e. xo(to) = XN(to), yo{to) = Z/iv(^o)- Assume further that {xN,yN) is a 
simple root of the system of equations Re{ip) = Im{ip) = 0. Eq.([T]) implies that ipoito) = 0, 
but not all the coefficients aio, Oqi, &io, &oi vanish at t = to- A third requirement is that 
the field current j = [Re{ip)'V{Im{ilj)) — Im{ip)'V{Re{ilj))]/{2i) should be divergence- free, 
V ■ j = 0, at the position of the nodal point. This follows from the continuity equation 
dp/dt + V ■ j = 0, since any zero of the wavefunction = is a local spatio-temporal 
minimum of p = l^/'p, thus dp/dt = at (0:0(^0)5 yo(^o))- The condition V ■ j = implies 
that 

^02 = "~<^20) ^02 = ~&20 ■ (3) 



Substituting the expansion ([T]) into ([2]), for t = to, and taking into account the previous 
conditions, the equations of motion take the form (up to second degree): 

1 



du 



G 



X 



(floi&io - aiQhQi)v + ^(ao2&io - aio&02)M^ + (^ao2&io - ^010602 - 011601)^'^ 



+ (ao2&oi - 001602)^ V + . 
dv 1 



(4) 



X 



aio^oi - ambiQjU + ^(aoi&o2 - 002^01)^^ + (^aoi&02 - ^ao2&oi - aii&io)^^ 



+ (aio&02 - ao2&io)M w + 



-K 



with 



G = (a^o + 6io)m^ + (floi + ^01)^^ + 2(aoiaio + &oi&io)^i + ■ ■ ■ (5) 

Equations (jlj) yield the ensemble of instantaneous flow lines (phase portrait) in the se- 
lected frame of reference for t = t^. Two questions are now examined, namely a) the typical 
form of the instantaneous phase portrait, and b) whether adiabatic conditions are satisfied, 
ensuring that the form of the phase portrait changes in time slowly, relative to the typical 
velocities along the particles' trajectories within this portrait. 



B. Phase portrait: nodal point - X-point complex 

The adiabatic approximation for t = to holds in space domains in which du/dt, dv/dt are 
large compared to the time derivatives of the coefficients aij, bij, and of the velocities Vx,Vy. 
We can then 'freeze' Oij, bij, Vx, Vy to their fixed values at t = to and treat the fiow (jl]) as 
autonomous. In such an approximation we find the following features of the instantaneous 
phase portrait around the nodal point {u,v) = (0,0): 



1. Nodal point 

In polar coordinates u = Rcos(j), v = Rsm(j) Eqs.(j4l) take the form: 

dR C2R^ + csR^ + caR^ + ... dcj) do + diR + diR^ + ... 



(6) 



dt G dt G 

where the coefficients Cj and dj depend on a) the coefficients a^j, b^j, b) the velocities (Vx, Vy), 
and c) powers of the trigonometric functions sin0, cos0. The lowest order terms of G, given 
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by Eq.([5]), read: 



G = R' 



^10 + ^lo) cos + (floi + &oi) sin + (o-oiciio + ^oi^io) sin: 



+ 0(i?3) . (7) 



The quantity in the square bracket of ([7]) is always positive. Thus, the second of Eqs.(l6]) 
imphes that for R small has a sign independent of 0, namely the same as the sign of the 
coefficient do = aio^oi ~ c^oi^io- Generically we have do 0. This implies that describes 
rotations clockwise, if do < 0, or anticlockwise, if do > 0. Furthermore, the angular frequency 
near the nodal point scales as = 0{R~'^). 

The flow lines of Eqs([6]) close to the nodal point (for R small) are given by 

d(t)~ do + diR + d2R'^ + ... ~ do \do dl ) "' ^' 
The coefficient C2 contains only terms of third degree in the trigonometric functions sin0, 
COS0. Thus, averaging Eq.Q over periods of the angle (which is fast for R small, = 
0{l/R^)) yields 

^=</3>i?^ + ... (9) 



02^1 

"o 



where the coefficient < /a > is given by 

< /a > (ay, K, Vy) = ^ ( ^ 

with i + j = 0, 1, 2. As explicitly demonstrated in Appendix A, for any non-zero value of the 
velocity of the frame of reference (\4, Vy) ^ (0, 0) we have < /a 0. Then the solutions 
of Eq.([9]) define spirals, i.e. 

^^^^^ T^2i?2</3>(0_0j • 



This is a spiral terminating at i? = 0, i.e., at the nodal point, when ^ oo (if < /s > < 0), 
or ^ — oo (if < /s > > 0). Thus, depending on the sign of < /a > the nodal point 
is either an attractor or a repellor. The only exception is when {Vx,Vy) = (0,0), i.e. the 
motions are considered in the rest frame. In that case we have < fj >= for all j > 3, 
i.e. the nodal point is a center. This follows trivially from the condition V ■ j = implying 
that if we set = — j j^du, the components of the current are given by = dH/dv, 
= —dH/du, which is equivalent to a Hamiltonian system {du / ds , dv / ds) = {ju,jv) under 
the non-uniform time parametrization ds = G^^dt = \ilj~'^\dt. Thus, in the rest frame the 
nodal point cannot be the limit of a spiral but it is a center of the instantaneous flow (see 
also [47,]). 



2. X-point 



The X-point {ux, vx) is a second critical point of the instantaneous flow found by setting 
du/dt = dv/dt = at {u,v) = {ux,Vx)- We flnd: 

]4 ^ Avx + Biuj^ + Civl + Diux vx + ■■■ 

Vy ~ -Aux + B2U\ + C2V\ + D2UX vx + ■■■ 

where the coefficients A, _Bj,Cj and Di are readily derived from Eqs.(jl]). If [ux^vx) are 
small, we flnd an approximate expression: 

Vx Vx 



Vy UX 

which, upon substitution to the flrst of Eqs.(j4]) yields a second order equation for, say, ux- 
The non-zero solution reads: 

giVxVy _yx , . 

''''' 92V^ + g,Vy^ + g,VxVy + g,V.! + g,V^Vy + g,VxV.f ^ \//'' ^ ' 
where all the coefficients gi depend only on the coefficients ajj, hij. In particular: 

Qi = 2(aio6oi ~ c^oi^io)) g2 = 2aii6oi + ^10602 ~ '^02^10 ~ 2aoi6ii, gs = 010^02 ~ 

5(4 = 2ao2&oi - 2aoi6o2, g5 = 2{al^ + bl-^), fi-e = -4(aoiaio + &01&10), gj = 2{alQ + bl^) . 

In numerical applications, Eq. (|T3l) is used to flnd a good initial guess for the position of 
the X-point, while better approximations are found by successive iterations of a root-flnding 
algorithm (e.g. Newton- Raphson) for the roots of the system of equations ([2]). 

The linearized equations of motion around {ux,vx) correspond to the linear system 
formed by the Jacobian matrix J[{dS/du, dS/dv), {u, v)] which is a 2 x 2 symmetric matrix 
with constant coefficients. Thus the eigenvalues Ai, A2 are real. In the limit of small ux, vx, 
the eigenvalues have opposite sign, since one readily finds that A1A2 = —A^/Gl + 0{1/R^), 
where G2 = 0{R^) is the quadratic part of G (Eq.Q). Hence {ux,vx) is an X-point with 
one unstable and one stable direction. Finally, the measure of both eigenvalues scales as an 



inverse power of the distance Rx = \/u'j^ + vj.. Numerically we flnd (see EKC) that the 
power-law scaling is A ~ Rx^ with p ~ 1.5. 



The numerator of the flrst of Eqs.( fT3i) is a O(V^) quantity, V = \/V^~-\-V^, while the 
denominator contains both 0{V'^) and 0{V^) quantities. If V is small we have Rx = 0(1). 
This, as shown below, sets an upper limit of validity of the adiabatic approximation for Rx- 
On the other hand, if V is large we have Rx = 0{V~^). 
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C. Conditions of validity of the adiabatic approximation 



Conditions for the validity of the adiabatic approximation are now visuahzed with the 
help of Figure 1 (schematic), showing the nodal point - X-point complex as viewed in a frame 
of reference of arbitrary velocity {Vx(t), Vy(t)) at two nearby times t = to (solid lines) and 
t' = to+At (dashed lines). Rx is the distance of the X-point from the nodal point, while ARq 
is the distance traveled by the nodal point with respect to this particular frame of reference 
within the time interval At. The vector v refers to the velocities of the particles' orbits 
as seen in the moving frame of reference. In the adiabatic approximation the velocities v 
must be large enough so that the flow integral curves change slowly relatively to the change 
of a particle's position within the time interval At. Since G in Eqs.Q depends, to the 
lowest order, quadratically on u, v (Eq.Q), the particles' velocities become arbitrarily large 
if u,v become sufficiently small. The linear velocities at a distance R from (xo,|/o) scale 
as |f| ~ R~^. Very close to (xo,yo) the motions are spiral-like with a frequency of order 
uj ~ \v\/ R ~ R~^i or period T ~ R^. The linear size of the 'nodal point - X-point complex' 
can be estimated by the distance Rx which is of order Rx ~ V^^ (Eq. f[T5]) ). The shift of 
the nodal point ARq within one period is estimated as ARq ~ T|y — Vo|, where Vq is the 
velocity of the nodal point in the rest frame. For the adiabatic approximation to hold, the 
following two conditions are sufficient and necessary: 

a) The shift ARq must be small with respect to the linear size Rx of the nodal point - 
X-point complex. This condition yields A_Ro << Rx, or -Rx|V — Vo| << Rx implying 



W -Vf 



^«1 . (14) 



V 

Thus the first condition is that V ~ Vq, i.e. the frame velocity V should be close to the 
velocity Vq of the nodal point with respect to the rest frame. 

b) The characteristic velocities within the 'nodal point - X-point complex' (i.e. for R < 
Rx) must be much larger than the rate of change of the coefficients aij,bij. Generically, 
the rates of change of aij,bij are in general 0(1) quantities (depending on trigonometric 
functions of the time and on the wavefunction's normalized amplitudes, see for example 
EKC). Thus \v\ should satisfy l-ul > 1, or, since \v\ ~ R^^ > Rx^, 

Rx ^V^^ <1 . (15) 
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FIG. 1: Schematic representation of the 'nodal point - X-point' complex as viewed in a frame of 
reference of arbitrary velocity {Vx{t), Vy{t)) at two nearby times t = to (solid lines) and t' = to + At 
(dashed lines). Ai?o is the length traveled by the nodal point within the time step At. Rx is the 
distance from the nodal point to the X-point at t = to- 

Thus, the second condition is that the velocity of the moving frame of reference should be 
large with respect to the rest frame. In particular, the rest frame itself represents a frame 
in which the use of the adiabatic approximation is not, in general, valid. 

D. Local Lyapunov exponents of scattered orbits 

The close encounters of the orbits with the nodal point - X-point complex can be approx- 
imated as scattering events, in which an orbit approaches the complex in a direction close to 
the stable manifold of the X-point and recedes from the complex in a direction close to the 
unstable manifold of the X-point. A simple model to describe the scattering is obtained by 
noting that if (without loss of generality) the axes are rotated so that at t = to the velocity 
of the nodal point is along the x-axis, i.e. xq ^ 0, i/o = 0, an expansion of the form ([T]) in 
the new coordinates yields the equations of motion 

du —Av + . . . . dv Au + . . . 

IE ~ Bu^ + 2Cuv + Dv^ + ... ~ H ~ Bu^ + 2Cuv + Dv^ + . . . ^ ' 

where A = aio&oi ~ Oqi&io, B = o^q + fe^O' ^ ~ '^oi'^io + ^oi^ioi D = al^+ b^^. The quadratic 
form in the denominator is always positive definite. The X-point is on the f — axis, with 
{ux,Vx) = {0, —AD^^Xq^), consistent with the scaling ~ = \xo\~^ found in the 
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previous subsection. Ignoring the higher order terms in the numerators of Eqs. (fT6|) causes 
the nodal point to become a center rather than the hmit of a spiral. This poses however no 
problem to the study of orbits being scattered by the nodal point - X-point complex since 
these orbits avoid penetrating the interior of the separatrix domain, close to the nodal point 
(see next subsection). Eq. lITBi) is then suggestive of the following simple model 

du —v.dv u . , 

to describe the scattering of orbits by the nodal point - X-point complex. 
The flow (ITTIl admits the integral 

which is obeyed by the scattered orbits locally, as long as the scattering lasts. The time 



evolution of i? = + t>2 is given by RdR/dt = XqB? — ihiR — \ InC)^. Far from the 
complex Eqs. llTTl) take the form du/dt ^ — xo, dv/dt ~ 0. Thus d{B?)/dt ~ —2xq{uq — xqI), 
where Uq = Rq is the initial condition of an orbit on the u-axis {vq = 0), implying 

^ Rl- 2xoUot + xlf . (19) 

The orbits passing outside the separatrix loop of the nodal point - X-point complex can be 
divided into 'type I' and 'type IF orbits (Figure 2a). Type I orbits surround the separatrix 
loop, while Type II orbits pass below the X-point, not surrounding the separatrix loop. In 
either case, the average time of a scattering event can be estimated by the non-trivial root 
for t of Eq. fll9l) with R = Rq, namely 

tscatter = O ( — — ] . (20) 



Deviations of t scatter from the estimate of Eq. fl20l) take place when the orbits are very close 
to the invariant manifolds of the X-point, since tscatter ^ C)0 as the initial conditions tend 
to a point on the stable manifold S. As shown in the Appendix B, such deviations lead to 
increased local Lyapunov characteristic numbers of the scattered orbits. The time evolution 
of the length of the deviation vector ^{t) = (Au, Av) has a characteristic 'profile' along the 
scattering, different for type I or type II orbits (Figure 2b). In the case of type I orbits, 
^{t) exhibits a rise and fall phase during the description of the separatrix loop. This phase 
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(a) (b) 




FIG. 2: (a) The thick looped curve is the separatrix passing from the X-point of the model (jl7p 
when xq = 3. The stable and unstable manifolds are marked S and U respectively. The stable 
manifold crosses the line u = 1 at Vg = —0.7785019.... The upper and lower thin solid lines represent 
a 'type 1' orbit (initial conditions u = l,v = Vg + 0.01) and a 'type 11' orbit (initial conditions 
u = l,v = Vg — 0.005) respectively, (b) The time growth of the deviations ^(t) for the type I (thick 
curve) and type II (thin curve) orbits. In both cases the initial deviation vector is taken = (1) 0). 
The rightmost vertical dashed line at t = 1 corresponds approximately to the time when the orbits 
reach u = —1, i.e. a position symmetric to their initial conditions with respect to the v = axis. 
The two left vertical dashed lines mark the time window 0.45 < t < 0.55 within which the type I 
orbit forms part of a loop around the nodal point. 

lasts for a time tioop estimated as tioop ~ -R^ ~ |io|~^- If l^^ol >> 1 we have tioop << t scatter- 
Most of the growth of C, takes place after t = Uoop, as the orbit recedes along the X-point's 
unstable manifold. In the case of type II orbits the ^{t) time profile exhibits a continuous 
rise from the start and ^(t) tends to stabilize after t = tgcatter- 

A theoretical quantitative estimate of the local value of the Lyapunov characteristic 
number in a scattering event is made in Appendix B. The growth of deviations is modeled 
by calculating the differential velocity of motions in two nearby integral curves of the flow 
( IT71) close to the X-point's stable and unstable manifolds. This modeling yields the scaling 
law 

f ~ (21) 
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FIG. 3: (a) The final value of (,/£,o (on a logarithmic scale) as a function of the initial condition 
V = vi of the orbits of the model ()17p taken on the line u = 1 at t = 0. The left and right curves 
correspond to xq = 3 and xq = 30 respectively. The vertical dashed lines mark the position v = Vg 
at which the X-point's stable manifold S crosses the line u = 1 in each case, (b) (,/£,o as a function 
of 6vi = \vi — Vs\ on a logarithmic scale. The two top curves correspond to xq = 3 (upper curve for 
type II orbits, lower curve for type I orbits). The straight solid lines passing through these curves 
represent a power-law fitting (,/$,o = A6v^^ with b = 0.95.. for the upper curve and b = 1.01 for 
the lower curve. The bottom two curves correspond to xq = 30. (c) A power-law fitting of the 
constant A as a function of xq {A = 0.89xq "'^'^ and A = 2.33i;Q "^'^^ for the lower (type I) and upper 
(type II) curves respectively). 

where and are the lengths of the deviation vectors of a scattered orbit before and after 
the scattering respectively, and 6vi is the initial distance of the orbit from the X-point's 
stable manifold far from the complex. The latter quantity is called the 'impact parameter'. 
The theoretical prediction (121 p is well reproduced numerically, by taking many trajectories 
in the model (fT7|) . for different values of xq (Figure 3). 

Since the local eigenvalue A of the linearized flow near the X-point scales as a positive 
power of Xq, Eq. (l2ip implies that the chaotic scattering takes place mainly in encounters 
in which A is relatively small (though non zero). This appears at first counter-intuitive. 
However, even if a trajectory is started close to the asymptotic manifolds of the X-point, it 
is in general far from the X-point itself, except for a short transit time of order ~ l/^o- Thus, 
to describe the total scattering correctly one has to take into account nonlinear terms of 
the expansion of the equations of motion, which introduce large deviations from the locally 
hyperbolic dynamics. On the other hand, the whole previous analysis relies on the use of the 
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adiabatic approximation, which, according to subsection II B, holds better when xq is large. 
For a nodal point - X-point complex to cause effective chaotic scattering, we thus have both 
an upper and lower restriction to the values of Xq. Precise upper and lower limits on xq, 
or, equivalently, the size of a complex Rx ~ l/io; ioT the complexes to produce effective 
chaotic scattering, can only be found by numerical experiments, as substantiated by specific 
examples in subsection III D. 



III. NUMERICAL STUDY 



A. A numerical example of the nodal point - X-point complex 

In our previous paper (EKC) we studied the 'nodal point - X-point' complex in a system 
of two harmonic oscillators 



H = l{pl + pl) + l{x' + icyr) (22) 



when the guiding field is the superposition of the ground state and the two first excited 
states [281] 

i;{x, y, t) = e"^^-^^ (i + axe-'' + hc^'^ye-'^^^"^') (23) 

while the frequencies are incommensurate, c^i = 1, u;2 = c = \/2/2. If we select a moving 
frame of reference such that its center (xo(t), ?/o(^)) coincides at all times with the moving 
nodal point, Eqs.(jl]) take the form: 

du 6c^''^i; sin(l + c)t 



dt G 

dv hc^/'^u sin(l + c)t — ahc^/'^v? sin ct 



xo (24) 
2/0 



dt G 
where G = + + G4 with G2 = (u'^/xl) - 2bc^^'^uv cos{l + c)t + b'^cxlv'^, G3 = 
— (26c^/^/a;o)M^f cos(l + c)t + 26^cxoUf^, G4 = lP'cv?v^ and xo5?/o are found by differenti- 
ating xo(t),yo(t), which are given by 

, , sin(l + c)t , , asint 

Figure 4 shows examples of the form of the nodal point - X-point complex in the above 
system, in a specific time interval. In all cases four asymptotic manifolds start from the 
X-point along pairs of opposite, stable or unstable, directions. One asymptotic manifold 
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goes towards the nodal point in a spiral way and the other three extend to infinity. The 
nodal point itself acts as an attractor or a repellor, and an asymptotic curve starting at 
the nodal point forms a spiral outwards. The sense of motion around the nodal point is 
determined by the sign of < /a > (Eq.([n])), which is given in this case by ^ 

^ / 1 + b^cxl \ 
^ ^ (^46ci/2a;4sin(l + c)toy ^ 

In view of Eq.f llOl) . if < /s >< the value of R decreases towards R = with increasing 
0, therefore the nodal point is an attractor if the spiral is described counterclockwise and a 
repellor if it is described clockwise. The opposite is true if < /s > > 0. On the other hand, 
the coefficient do in Eq.([6]) turns out to be equal to do = sin(l + c)to (while G is positive for 
any {u,v) ^ (0,0). Thus, if sin(l + c)to > (or< 0) the spiral is described counterclockwise 
(or clockwise). When < /a >= (and sin(l + c)to 7^ 0) the nodal point changes from an 
attractor to a repellor. Then we have a Hopf bifurcation, followed by the formation of a 
limit cycle. 

As an example we consider the evolution of the manifolds between t = 1.25 and t = 1.35 
(Fig. 4). For t = 1.25 (Fig. 4a) we have < /s > < and sin(l + c)t > therefore the nodal 
point is an attractor and the orbits appoaching it are spirals described counterclockwise. 
One orbit of this type is one of the unstable manifolds of the X-point. The other unstable 
manifold and the two stable manifolds of the X-point extend to infinity. In particular, 
the upper stable manifold escapes downwards after making an almost complete rotation 
(backwards in time) clockwise around the nodal point. 

The attraction of the orbits by the nodal point terminates when a transition of < /a > 
takes place from negative to positive, near t = 1.294. Then, the nodal point becomes 
a repellor and a limit cycle is formed around it. The limit cycle moves outwards (e.g. 
t = 1.296, Fig. 4b). The orbits both inside the limit cycle (i.e. close to the nodal point) and 



^ In our previous paper (EKC) an error appears in the second factor of < /s > in Eq.(35). This equation 
is derived from Eq.(A4) which is correct. However, in Eq.(35) the term xqXq should be replaced by 
^"^(i+b^cJ^)"^ ■ numerical results change only slightly. Note that Fig. lib of EKC gives only the second 
factor of < /a >, that changes sign. Three more typos have been found in the caption of Fig. 12, namely 
the integer part of to in the cases (a),(b),(c) is 175, as in case (d). 
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FIG. 4: The form of the nodal point - X-point complex in the EKC model (equations of motion 
given by ()24p with a = b = 1, c = for the times indicated in panels (a) to (d). The nodal 

point is an attractor at t = 1.25 (a). A Hopf bifmxation takes place near t = 1.29415. The nodal 
point becomes a repellor and the associated limit cycle moves outwards at subsequent times. E.g. 
at t = 1.296 it has the position shown in (b). At t = 1.303 the limit cycle reaches the X-point 
(c). Then the limit cycle disappears and all the orbits inside the nodal point - X-point complex 
are repelled away from the complex (e.g. at t = 1.35,(d)) 

outside the limit cycle (between the limit cycle and the X-point) are attracted by the limit 
cycle. The orbits can enter the complex only via a very narrow channel formed by the two 
stable manifolds below the X-point. As the limit cycle moves outwards and approaches the 
X-point, any orbit approaching the limit cycle is dragged closer and closer to the X-point. 

The limit cycle reaches the X-point at t = 1.303 (Fig.4c). Then the unstable manifold 
from the right joins the upper stable manifold and together they form a separatrix. For still 
larger t (e.g. t = 1.35, Fig.4d) the limit cycle has disappeared and the upper stable manifold 
approaches the nodal point via spiral rotations (backwards in time). On the other hand, the 
two unstable manifolds go to infinity on the left, one directly, and the other after an almost 
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FIG. 5: The value of < /s > as a function of the time t in the interval < i < 10 in the EKC 
model with a = b = 1, c = \/2/2. The dots mark the points where < fs >= 0. 

complete rotation (counterclockwise) around the nodal point. 

The transition displayed in Fig.4 constitutes a Hopf bifurcation. Hopf bifurcations take 
place whenever < fs > crosses a zero value. The bifurcation displayed in Fig.4 is called 
direct (the limit cycle is formed first near the nodal point and later it disappears at a 
separatrix). However, inverse Hopf bifurcations are also commonly observed, in which the 
limit cycle moves towards the nodal point. The rate of appearance of direct or inverse Hopf 
bifurcations is a few per period (which is of order 27r). A typical survival time for limit 
cycles is AT ^ 10"^ (in Fig.4 we have AT = 0.008). 

The value of < /s > follows a time evolution as exemphfied in Fig. 5. The value of 

< /a > becomes infinite when sin(l + c)to — and when sinc^o = 0, with to ^ 0. In Fig. 5 

< /a >= oo at times to = kn/{\ + c) with k = 1 . . .5 and to = tt/c and 27r/c. Between 
k = 1 and 2 the sign of sin(l + c)to is negative, thus the nodal point is a repellor whenever 

< /a > > 0. Between k = 2 and 3 the sign of sin(l + c)to is positive, thus the nodal point 
is a repellor whenever < /s > < 0, and so on. Between two points with < /a >= oo there 
may be a number of times where < >— (one, two, or three times in Fig. 5). Evolutions 
of the phase portraits similar to Fig.4 appear close to all the time moments when < fs >— 0. 

Using the above rules, in the time interval < t < 10 the nodal point is found to be a 
repellor for about 47% of the total time. Whenever the nodal point is a repellor it cannot 
in general be approached by any orbits. A rare exception is the case in which an inverse 
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Hopf bifurcation takes place at the nodal point. Then the orbits which are initially in an 
extremely narrow channel of the flow, formed between the two stable manifolds outside the 
X-point, approach the limit cycle which tends to the nodal point. Such events can only last 
for times AT ^ IQ-^. 

Similarly, when the nodal point becomes an attractor, only an extremely narrow channel 
of the flow formed by the stable manifolds allows for the orbits to go deeply inside the 
complex, i.e., to approach the nodal point. But this channel also disappears in transient 
time intervals AT ^ 10~^ in which the nodal point is protected by a limit cycle. 

We conclude that the asymptotic curves of the X-point in most cases do not allow close 
approaches to the nodal point. Only if two conditions are satisfied, namely that (a) the 
inner asymptotic curve is unstable, and (b) the nodal point (or a limit cycle approaching it) 
is an attractor, we may have close approaches to the nodal point. In all cases, however, the 
approach is only possible for a set of initial conditions of extremely small measure, i.e., the 
orbits in general avoid the nodal point. 



B. Type I and Type II interactions 

Figures 6 and 7 show now examples of Type I (Fig. 6) and Type II (Fig. 7) interactions 
of a quantum trajectory with the nodal point - X-point complex in the above system. The 
trajectories are viewed in the moving frame of reference {u, v) centered at the nodal point and 
they are superposed to the background instantaneous velocity flow at the indicated times 
t. The deviation vector ^ is calculated numerically by the variational equations of motion. 
Its local direction is indicated by the thick arrows in each panel, while the time evolution 
of the normalized length ,^/,^o is shown in the last panels of Figs. 6 and 7. These two curves 
are compared to the theoretical curves of Fig.2b for type I and type II events respectively. 
Note in particular that in the case of the type I event the growth of ^ during the description 
of the first half-loop is nearly compensated by a decrease in the second half-loop, and the 
deviations start growing essentially after the description of the loop. The two local peaks of 
^{t) in the interval of the loop description can be explained by the fact that the loop is not 
perfectly circular (see [48] for an explanation of the behavior of ^ in non-circular invariant 
curves). Also, in the case of the type II orbit the growth of ^ does not take place very close 
to to the X-point, but all along the scattering event, in accordance with the theory of section 
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FIG. 6: A type I event in the EKC model with a = 6 = 1, c = \/2/2, and initial conditions 
of the trajectory xq = 1.71510, yo = 1.29285 at t = 0. In panels (a) to (e), the nodal point 
appears as a gray thick dot at (0,0), while the X-point appears as a black thick dot (denoted X). 
The background small arrows indicate the local instantaneous velocity flow in the {u, v) frame of 
reference. The thick long arrow shows the direction of the deviation vector ^ for the same orbit 
(initial conditions = (1)0))- Panel (f) shows the time evolution of the normalized length ^/^o 
on a logarithmic scale. 



II. 



C. Multiple nodal points 

In EKC we considered the orbits in simple ■?/'— fields of the system of two harmonic 
oscillators given by the Hamiltonian fl22|) . The eigenfuctions are 



where the energy of the state {nin2) is 



En-,n2 = ( ^ + '^l ) ^1 + + ^2 ) (28) 
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FIG. 7: Same as in Fig. 6 for a type II event of the same orbit. 

and Hn are Hermite polynomials. The following eigenfunctions are explicitly referred to in 
the sequel: 

^00 = e-^^^+'^^^e ^ , = e's^^+'^^'e ^ x, = e'^^^+'^^^e ^ x^fcy 
vl/20 = e-i^'+'^'e~^-^ {x^ - l) , ^>,o = e-^^^+'^^^e"^^ (x^ - 3a;) . (29) 

In EKC we considered particular quantum trajectories in the case ^ = \E'oo + '^^oi + ^c^/^\E'ii 
(Eq. (l23l) ). with a,b real, yielding one moving nodal point. Here we shall consider the cases 
with a progressively higher number of moving nodal points in the same Hamiltonian model. 
We call 'nodal lines' the trajectories of the nodal points. For the nodal points to be moving, 
there are restrictions on the choice of combination of the eigenfunctions, since for particular 
combinations there are no nodal lines but isolated nodal points appearing at specific times 
only. For example, if the wavefunction consists of the sum of three eigenfunctions = 
^nin2 + '^^n'n', + n" n" ^ith equal quautum numbers 77.2 = rig = n'2, the nodal points 
\E' = satisfy the equations 



^ = H^, (x) + ae-*V"i-"^i*i/„, (x) + fee'^"^ ""'^'H^" (x) = . (30) 



20 




FIG. 8: Nodal lines of the wavefunction ^ = *oo + «*io + ^^ii + e^20 when a = 6 = 1, c = \/2/2, 
e = 0.1. 

Thus we have two distinct equations for the real and imaginary parts of implying that 
we have solutions for x only at specific times t. The same happens if ni = n'^ = n[. The 
same is true if we have more terms in but with the same quantum number ni, or ?t,2, in 
all the terms. Such cases are not examined below. 

The nodal lines may enclose domains in the configuration space devoid of nodal points 
for all times. In such domains the quantum trajectories turn to be regular. Such empty 
domains are found even if we take combinations of eigenf unctions yielding simultaneously 
more than one nodal points. Examples are: 

1. Case = ^'oo + a*io + + e^20 

If we add a fourth term of the form \l/2o in Eq. (!23!l . i.e.: 

= ^00 + a^io + b^u + e^2o (31) 
the real and imaginary parts of the equation \1/ = take the form 

1 + ax cos t + hx\/cy cos (1 + c) t + e [x^ — l) cos 2t = 

ax sm t + 6a;yc?/sin(l + c)t + e(x^-l)sin2t = . (32) 
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Multiplying the first equation by sin(l + c)t and the second equation by cos(l + c)t and 
subtracting, we find the equation 



sm 



'1 + c)t + ax sin ct + e (x^ — l) sin (c — 1) t = 



with solution 
1 



X 



(33) 



1 /2 

-—{—a sin ct ± Ta^ sin^ ct — 4e sin (c — l)t (sin (1 + c) t — e sin (c — 1) t)] } 
2e sin (c — 1) t 

(34) 



If sin ct is not close to zero and e is small we find 

1 



X 



-{—a sin ct ± a sin ct 



+ 



2e sin {c — l)t 
2e2sin2(c- l)t 2e2sin2(c- l)t sin^(l + c)t 



^ 2e sin(c — l)i(: sin(l + c)t 
sin^ ct 



} 



(35) 



sin^ ct sin^ ct 

If xo denotes the solution ( l25i) . for e = 0, the solution close to x = xq for e small is the one 
with the plus sign 



sin(l + c)t esin(c— l)t 

X = : : h 



a sm ct 



and 



a sin ct 
1 



sin^(l + c)t 
sin^ ct 



e sin(c — l)t ^ 
a sin ct 



1 — X, 
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[— axsint + e(l — x^) sin2t] . 



(36) 



(37) 



x6ycsin(l + c)t 

When t = fcvr we have y = and x = ±^ + O(e^). In particular if a = 1, we have x = ±1 
exactly. This is seen in Fig. 8. Besides this solution we have y = also if —ax + 2e(l — 
x^) cost = 0, hence 

a I 16e^cos^t 

X = — 

4e cos t 

and if we take terms up to 0(e) we find 

2e 



1± A 1 



X = — COS t and x 
a 



2e cost 



(38) 

2e cos t a 

These solutions must match the solution (!36|) and this matching should give the time t. The 
first solution of fl38l) is of 0(e) and cannot ever match the solution fl36l) or flMl) which is of 
0(1). On the other hand the second solution of (155]) is of 0(l/e), i.e. a large number, and 
this can match the solution fl5B]) if Xq is large. In fact in Fig.8 we see that there are solutions 
with y = for |x| > 8. 

In conclusion, similarly to the case considered in EKC, for e small the nodal lines leave 
an empty central domain, delimited by hyperbola-like boundaries. 
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FIG. 9: Nodal lines of the wavefunction ^' = ^'oo + a^'2o + t^ii when c = ^2/2 and (a) a = 1.23, 
b= 1.15, (b) a = 1, 6 = 1.15. 
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FIG. 10: (a) Nodal lines of the wavefunction ^ = '^qq + 0^*30 + when a = 1.23, h = 1.15, and 
c = \/2/2. (b) The time evolution of XQ{t) for either one or the three solutions of 

2. Case = ^'oo + «*20 + ^^ii 

In the same way as above, the real and imaginary parts of \I/ = yield 



1 + a [x^ — 1) cos 2t + bx\^y cos (1 -|- c) t = 
a (x^ — 1) sin 2t + hx\fcy sin (1 -|- c) t = . 



(39) 



Multiplying the second equation by cot (1 -|- c) t and subtracting from the first we find 

sin (1 + c)t + a (x^ - 1) sin (c — 1) if: = (40) 
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hence 



X 



2 



a [X 



2-l)sin2t sin2t 



6sin(c+l)t 6sin(c— l)t 



(41) 



a sin (c — 1) t 

Equations (14T|) give two nodal points with opposite x and y whenever the restriction x^ > 
is satisfied, and no nodal points when it is not. Thus, the nodal points exist only in particular 
time intervals. 

In general the nodal lines enclose an empty region near the origin (Fig.9a). However in 
the particular case a = 1 the nodal lines reach the center x = y = (Fig. 9b). In fact, y = ii 
sin2t = 0, i.e. t = 0, n/2, n... The corresponding values of are : For t = 0, = 1 — ^^^^-^ , 
for t = 7r/2, x^ = 1 + ^, for t = tt, = 1 — ^ etc. The last solution exists only if a > 1, and 
it yields x = if a = 1 exactly. On the other hand the solution x^ = 1 + ^ always exists if 
a > 0, meaning that the nodal lines intersect the x— axis at the points x = ±(1 + 1/a)^/^ 
(Figs.9a,b). 



3. Case = ^'oo + a^'so + ^^ii 

In this case the real and imaginary parts of \E' = yield 

1 + a (x^ — 3x) cos 3t + hx\fcy cos (1 + c) t = 

a {x? - 3x) sin 3t + hx^fcy sin (1 + c) t = . (42) 

Multiplying the first equation by sin (1 + c) t and the second equation by cos (1 + c) t and 
subtracting we find 

sin (1 + c) t + a (x^ - 3x) sin (c - 2) t = (43) 

hence 

o sin (1 + c) t sin 3t , , 

x^ - 3x + '— = 0, y = — (44 

a sm (c — 2) t Oa/cx sm (c — 2) t 

The third degree equation has three real roots if 



sin^ (1 + c) t 



< 1, (45) 



4a2 sin^ (c - 2) t 

otherwise only one root is real. 

The nodal lines in this case have the form of Fig. 10a, leaving again empty central domains 
in the configuration space. We have three nodal points if the inequality fHSjl is satisfied, 
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and one nodal point otherwise. In fact, there are distinct time intervals within which one 
of the nodal points comes from or goes to infinity (Fig. 10b). For example, after the time 
ti = 2.43, one nodal point (point 1) starts approaching the central region from x oo. A 
little later (at t = thij = 2.75), a pair of nodal points (2 and 3) emerge at a; ~ —1. At the 
time if: = t'^^^j = 4.54, point 2 joins point 1 nearly at x = 1, and after t = t'^^^j these two points 
disappear, while point 3 tends to x — — oo at t = ^2 = 4.85. Similar phenomena take place 
at subsequent intervals of time. The times etc. are called 'bifurcation times'. Such 

bifurcations are important for the level of chaos of the trajectories approaching the nodal 
points, because close to a bifurcation time the speed of the bifurcating nodal points (e.g. 
Xq) which enters into the estimates of local Lyapunov characteristic numbers (Eq. fl2ip ). is 
large (see numerical simulations below). 



D. The degree of chaos for ensembles of chaotic trajectories 

As a first example, we consider orbits in the ■j/'-field = \E'oo + ct^2o + b^^ii (subsection 
III C 2) when a = 1.23, b = 1.15, c = \/2/2. The nodal points in this field appear in pairs, 
within some time intervals. The X-points are calculated by a Newton-Raphson method with 
a precision tolerance of 10~^^, loading (fT3l) as initial guess values of {ux,vx)- Fixing the 
frame of reference on one of the nodal points, two X-points are found numerically. One 
X-point is close to the considered nodal point and the other is far from it. Nevertheless, 
the distant X-point is irrelevant to the dynamics, because in that case we have (for, say, the 
nodal point (xoi,?/oi)) \V ~ K)2| = \Voi - K)2| = 2|Koi| = 2|t^| (since V02 = -Vqi), implying 
1^ — Vo2|/|V^| = 2, i.e., a violation of the condition (IMI) . Thus, in the numerical calculations 
we only take into account approaches to the X-points found in the vicinity of each nodal 
point in its own frame of reference. 

For any trajectory with initial conditions x{0),y{0) there is a centrally symmetric orbit 
with initial condition — x(0), —y{0). Figure 11 shows three orbits: regular (Fig. 11a), weakly 
chaotic (Fig.llb), and strongly chaotic (Fig.llc). The degree of chaos is measured by the 
quantity 

(46) 

where ^{t) is the length, at time t, of a deviation ^ = (Ax, Ay) from an orbit j x{t),y{t) J , 
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FIG. 11: Examples of quantum trajectories in the tp— field * = ^oo + a^20 + b'^ii when a = 1.23, 
b = 1.15, and c = \/2/2. (a) An ordered orbit (initial conditions x(0) = —1.5, y{0) = 0.1275, (b) a 
weakly chaotic orbit (initial conditions x{0) = 0.850901842117, y{0) = 1.191571712494, and (c) a 
strongly chaotic orbit (initial conditions x(0) = -1.231356695294, y(0) = 0.840584903955. In all 
three panels the orbits are plotted up to t = 1000. (d) Time evolution of the 'finite time Lyapunov 
characteristic number' x{t) for the three orbits. 

calculated by the variational equations of motion. This is called 'finite time Lyapunov 
characteristic number' and the limit \imt^Qox{t) yields the Lyapunov characteristic number 
of an orbit. In the case of the regular orbit, the quantity x(^) (Fig-Hd) decreases as a 
power law x(i) ~ t~^. In the case of the strongly chaotic orbit, after some transient time 
the quantity x{t) decreases slowly and it tends to stabihze to a value x{t) — 5 x 10~^ at 
t — 10^. On the other hand, in the case of the weakly chaotic orbit, there is a temporary 
stabihzation of x(^) up to f = 10^, followed, however, by a decrease up to i = 7 x 10^. 
Beyond this time x(i) increases again up to the value x = 10~^ at t = 10^, showing no signs 
of stabilization. 
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In general, weakly chaotic orbits exhibit long transient intervals in which x fluctuates 
around values typically one order of magnitude smaller than the stabilization value of the 
strongly chaotic orbits. A careful inspection shows that in most cases this behavior of the 
weakly chaotic orbits can be characterized as a 'stickiness' phenomenon (see e.g. j^ol), 
namely the orbits behave essentially as regular in a transient time interval. In Fig. lib this 
tendency is observed for the weakly chaotic orbit, which, besides the chaotic oscillations, 
shows a domain of enhanced density similar to the domain filled by the regular orbit of 
Fig. 11a. 

Also, the difference in the evolution of for the three orbits is related to their frequency 
of encounters with nodal point - X-point complexes. The background points in Figs.lla,b,c 
show the distribution of the X-points in the configuration space, which remains practically 
unaltered after a time t = 1000. The X-points occupy domains similar to those occupied 
by the nodal points (Fig. 9a). Setting an upper threshold distance dmax, and splitting the 
time evolution of the orbits into time segments of width At, we may count the number of 
time windows within which a trajectory approached the X-point at a minimum distance 
d < dmax- In the numerical calculations we set d^ax = 0.2, and At = 0.1, and find a 
number of approaches, up to a time t = 10^, equal to N{d < 0.2) = for the regular orbit, 
N{d < 0.2) = 106 for the weakly chaotic orbit and N{d < 0.2) = 389 for the strongly chaotic 
orbit. One can check that the ratios of the values of N{d < dmax) of the three orbits remain 
practically invariant if another choice of dmax is made, provided that dmax is bounded (e.g. 
dmax does not exceed unity). 

In the system = \I^oo + ct^so + (subsection III C 3, one or three nodal points) we 
also find regular, chaotic and weakly chaotic orbits. In the system of subsection III C 1 we 
find orbits similar to those of EKC as long as the parameter e is small. 

Figure 12 shows now the main result. Having run 500 orbits with initial conditions 
taken randomly in the box — 1.5<x< 1-5, < y < 1.5 in the three systems (i) \1/ = 
\E'oo + ct^io + ^^ii (EKC, one nodal point), (ii) \1' = \I^oo + '^^2o + &^ii (two nodal points) and 
(iii) \E' = \E'oo + 0^30 + (three or one nodal points). Fig. 12a compares the distributions 
of the finite time Lyapunov characteristic numbers x{t = 10^) for the three ensembles of 
orbits. The time t = 10^ is long enough to extinguish transient effects of x{t) for uiost 
orbits. It is immediately clear that model (ii) yields a significantly larger degree of chaos 
than models (i) and (iii). Model (i) yields a bimodal distribution, with a large proportion 
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}(:(t=io') logio(Rx) 

FIG. 12: Histograms (line diagrams) of the number of orbits AN of which the finite time Lyapunov 
characteristic number at t = 10^ is in the interval [x — 0.0025, x + 0.0025), where x is the value 
shown in the abscissa. The orbits result from 500 initial conditions taken randomly in the box 
— 1.5 < X < 1.5, < y < 1.5. The solid line with squares corresponds to the EKC model (i) 
* = *oo + a^io + (one nodal point), the dashed line to the model (ii) = ^'oo + a*20 + 
(two nodal points), and the solid line with crossed circles to the model (iii) * = *oo + a*30 + 
(three or one nodal point). In all three cases a = 1.23, b = 1.15, c = V2/2. (b) Histo gram of the 
values of IoqiqRx (nodal point - X-point distance) when all the nodal point - X-point complexes 
are calculated at the times t = n x 10~^, n = 1, 10^. The dotted solid line refers to the EKC 
model (i), the dashed line to the model (ii), and the solid line to the model (iii), with parameters 
as in (a). 

of regular orbits (x = 0) and also a local maximum of the distribution of the chaotic orbits 
at X ~ 0.025. Model (iii), on the other hand, has a small number of perfectly regular orbits, 
but the main bulk of its chaotic orbits is also in rather small values of x (x < 0.01). 

In the systems (ii) and (iii) the double or triple nodal points appear only in certain time 
intervals. In a total time t — 10^, divided in segments At = 10~^, the total number of 
nodal point - X-point complexes detected in a box (x, y) e [—3, 3] x [—3, 3] are 8.4 x 10^ 
for the system (i), 5.1 x 10^ for the system (ii), and 1.47 x 10^ for the system (iii). Thus, 
chaos appears less pronounced precisely in the system exhibiting the largest number of nodal 
point - X-point complexes, i.e. the system (ii). This phenomenon can be understood if we 
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take into account the theory of subsection II D, and in particular the fact that the values 
of the local Lyapunov characteristic number have a 0{V~^) dependence on the speed of the 
nodal point (Eq. ipT]) ). or, equivalently, a 0{Rx) dependence on the nodal point - X-point 
distance. Plotting the histograms of the values of Rx for all three systems (Fig. 12b) renders 
immediately clear that in the system (ii) (two nodal points) the main bulk of the histogram 
is at values of Rx larger than in both the systems (i) and (iii) (one or three nodal points), 
i.e. the system (ii) has the more effective chaotic scatterers (complexes) from all three 
systems. In the system (ii), Rx is mainly distributed over the range 0.1 < Rx < 1 (with a 
mean < R^ >~ 0.2). Thus, both the condition of validity of the adiabatic approximation 
{Rx < 1, vertical dashed line in Fig. 12b) and the requirements for effective chaotic scattering 
{Rx large) are fulfilled. In the system (iii), the speed of bifurcating nodal points is large 
near the 'bifurcation times' tbf (see Fig. 10b and the relevant discussion), and this reduces 
the effectiveness of the chaotic scattering. In the case of the system (i) , there is a significant 
percentage of non-effective complexes {Rx < 10^^) or of complexes violating the condition 
of adiabaticity (i.e. with Rx > 1). In this system we thus find less chaos than in the system 
(ii), and also a large number of perfectly regular orbits. 

The value of x fo^^ particular orbits is in general an increasing function of the number 
of encounters with nodal point - X-point complexes, but this relation presents considerable 
scatter and also noticeable exceptions. In Figs. 13a, c,e the number of consecutive approaches 
of a trajectory to a nodal point X-point complex are estimated by the index N{d < 0.2) 
(number of approaches at a distance d < 0.2; similar results are found if N{d < 0.5) is used 
instead). The tendency of x to increase with is clear in all three systems. The exceptions 
refer to orbits in the lower right part of each of Figs.l3a,c,e. These exceptions disappear, 
however, if we use a corrected index for the number of encounters, by the requirement that an 
encounter is only counted provided that the size Rx of the complex during it is not very small. 
Figures 13b, d,f show the scaling of x versus the corrected index N{Rx < 0.5, d < Rx/2), 
for the same orbits, i.e. we count only the encounter events in which the distance Rx is 
not larger than twice the distance d at which the orbits have the closest local approach 
to a nodal point X-point complex (the factor two is rather arbitrary; in general we can 
set Rx > 0{d)). The second condition, Rx < 0.5, ensures that only complexes being 
well within the regime of validity of the adiabatic approximation are selected (this is also 
arbitrary; any limit Rx < 0{1) can be used). These extra conditions immediately yield a 
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FIG. 13: The finite time Lyapunov characteristic number x at t = 10^ as a function of the number 
N{d < 0.2) of approaches of an orbit to the X-point at a distance smaller or equal to 0.2 for 
the systems (i),(ii) and (iii), shown in (a), (c) and (e) respectively. For each system, 500 initial 
conditions are taken randomly in the box —1.5 < x < 1.5, < y < 1.5. The finite time Lyapunov 
characteristic number x versus the corrected index N{Rx < 0.5, d < i?x/2) (see text) is shown in 
(b), (d) and (f) respectively. 



lower number of recorded events for all the orbits than by the index N{d < 0.2). The main 
effect however is that all the exceptions of Figs.l3a,c,e disappear by utilizing the corrected 
index, thus yielding a better correlation of % with N{Rx < 0.5, d < Rx/2). 



IV. CONCLUSIONS 



We developed the general theory of motion in the vicinity of a moving 2D 'quantum 
vortex', i.e. a nodal point of the wavefunction, in the trajectory (Bohmian) approach of 
the quantum flow, and we discussed the origin and quantification of chaos for the Bohmian 
trajectories. Our main findings can be summarized as follows: 

1) The flow in the vicinity of a moving nodal point is non-autonomous, but under suitable 
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'adiabatic' conditions it can be treated as nearly autonomous. Two necessary and sufficient 
conditions are found: a) the equations of motion must be taken in a moving frame of reference 
centered at the nodal point, and b) the latter's velocity must be large in the rest frame. 

2) Developing an arbitrary wavefunction up to terms of second degree with respect 
to the distance from the nodal point, we demonstrate that the appearance of nodal point - 
X-point complexes is a generic feature of the configuration space. There are two stable and 
two unstable manifolds emanating from the X-point associated to each nodal point. One of 
these manifolds continues as a spiral approaching the nodal point, while two other form very 
narrow channels allowing communication with the interior of the complex. The nodal point 
undergoes consecutive Hopf bifurcations. Whenever a Hopf bifurcation takes place a limit 
cycle is formed around the nodal point for transient time intervals. As a consequence of all 
these facts, it is shown that most trajectories do not penetrate deeply into the complex. 

3) On the other hand, the chaotic orbits are scattered by the complex via encounters of 
'type r (forming a loop around the complex) or of 'type IV (no loop). A theoretical estimate 
is given of the local Lyapunov characteristic numbers in separate encounter events. The local 
Lyapunov characteristic number scales as an inverse power of the speed V of the nodal point 
and of the distance of the scattered trajectory from the X-point's stable manifold (impact 
parameter) far from the complex. The size of the complex (distance Rx of the X-point from 
the nodal point) scales as ~ V~^. The chaotic scattering is most effective when the speed 
of the nodal point is relatively small, or Rx is large. But Rx is also hmited by the extra 
condition of adiabaticity (Rx < 1)- Numerically, we find most effective chaotic scattering 
events taking place in the range 0.01 < Rx < 1, Rx ~ 0.1 being an optimal value. 

4) We provide numerical examples of the loci occupied by the nodal points and the X- 
points in different examples of superposition of a number of eigenstates in a 2D harmonic 
potential model. In particular, we examine three models with (i) one, (ii) two, or (iii) 
three nodal points, and identify the domains of each system devoid of nodal points. The 
trajectories having no overlap with the domains of nodal points turn to be regular. There 
are also weakly chaotic trajectories, exhibiting stickiness phenomena, and strongly chaotic 
orbits having a significant overlap with the domains of nodal points. The system with the 
smaller number of complexes (system (ii)) turns to have the largest degree of chaos. This 
is explained by examining carefully the properties of the complexes and demonstrating (on 
the basis of the theory of section II) that the most effective chaotic scattering events are 
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produced in the case of the system (ii). 

5) The 'finite time Lyapunov characteristic numbers x{t) of the trajectories have a nearly 
hnear correlation with the number of encounters with nodal point - X-point complexes, but 
with considerable scatter. The scatter is reduced, and most exceptions disappear when only 
'effective' events are counted. The effectiveness criterion takes into account the requirement 
that the trajectory approaches the X-point at moments when the complex is relatively large, 
implying that the chaotic scattering is strong. 
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APPENDIX A: EQUATIONS OF MOTION FOR A GENERIC tA-FIELD. 

The first of Eqs.Q is obtained by multiplying the first and the second of Eqs.dlj) by u and 
V respectively and adding the results. We then find u{du/dt)+v{dv/dt) = RdR/dt. Similarly, 
the second of Eqs.([HD is obtained by multiplying the first and the second of Eqs.Q by v 
and u respectively and subtracting the results. We then find v{du/dt) ~u{dv/dt) = Rd(f)/dt. 
After these operations, the values of the coefficients €2,03, do and di of Eq.([6]) are readily 
evaluated. Together with the conditions 092 = —0.20, &02 = — the average value of the 
coefficient /a: 

< /a > (flij, hj, Vy) = — — 

^TT Jq \do 

takes the form: 



C3 C2di 

d^ 
"0 



1 

< /s > 



( 2aoia^Q6o2+ctoi'^ii^io~'^io'^ii^io^2aioao2Ctoi^io+'^ii^oi^io 



4(aio&oi - aoi&io)^ 

+2aio6oi&02&io - 2ao2&oi&?o ^ '^n^io - floi'^io^H + aio^n " aio&oi^n + «io&?o^ii 



-Vy \^aoiao2aioboi — 2aioaoi^o2 + ^lo'^^ii^oi ~ 001*^11^01 + oii^io^oi 



+2ao2&io&oi ~ 2aoi6io&o2^oi ~ o-ubQi — ctioC^oi^ii + Oqi^ii " Ooi^io^n + '^oi^oi^ 
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+ (K - K/ ) ( «oi«io + «oi«io + «oi«io&oi + «oi^oi^io + a^oboiho + &oi^io + «oiaio&io + ^oi&io 

,4 „4 , o„2 t2 , t4 7,4 r,„2 7,2 



+V;V; Ooi - + 2aoifeoi + &oi - ^lo - 2aio&io 



Since all the terms in the above expression have either \4 or as a coefficient, it follows 
that < /a >= if Vx = Vy = 0, i.e. the nodal point is a center in the rest frame, and an 
attractor or repellor in any other moving frame of reference. 



APPENDIX B: LOCAL GROWTH OF THE DEVIATIONS ^{t) IN A TRAJEC- 
TORY - NODAL POINT - X-POINT SCATTERING EVENT 

Referring to the model f|T7|) of subsection II D, let T{C) be the time required for an orbit 
to traverse the complex along one of the integral curves given by (Eq lTSl) . For simplicity 
(and without loss of generality) we consider the case Xq > and identify T(C) to the time 
needed for an orbit starting on a curve C, at -u = 1, until the orbit crosses u = —1, namely: 

nC) = 2 / , (Bl) 



V 



where vo and vi correspond to the w— values satisfying 

i.e. the values of v on the curve C for m = and u = Ui = ±1 respectively. The X-point is 
located at 

Ux = 0, Vx = (B3) 

Xo 

and its asymptotic curves have the C-value 

C = Cx = 4^ . (B4) 

When Xo is large the asymptotic curves become nearly horizontal a little further from the 
X-point, i.e. |fi(C^)| is small. The same holds true for nearby curves with C ~ C^- The 
value of f 1 can then be found approximately by expanding ln(C) = 2xof i + ln(l + vf) to 
second order in vi, yielding 

MC)^- ■ (B5) 
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In view of (lB4p we have |lnC| ~ 21n|i;o| << l^ol- Thus, setting fi(C) ~ ln(C)/(2xo) is 
nearly always a sufficient approximation. 

Consider now the orbits on two neighboring integral curves C, C + AC such that T{C) > 
T{C + AC) = T(C) + AT. Far from the complex the velocity of the orbits is \V\ ~ —Xq. 
It follows that at the time t = T{C) the two orbits are a distance ~ |xoAT| apart. Thus 
the initial deviation dvi = has grown to ^ ^ + |a;oAT|. But AT = {dT/dC)AC and 
by virtue of flB5l) AC ~ 2Cxodvi = 2Cxo^o- Thus, the final value of the deviation can be 
estimated as: 

c) . (B6) 

Equation ( 1B6P states that the growth of deviations is proportional to the differential rate 
of description of the integral curves passing close to the X-point. If u > 1 the velocity 
stabilizes to ~ —Xq along all the integral curves. This explains the stabilization of ^ in 
Fig. 2b. Furthermore, \dT/dC\ increases as C tends to Cx, since limc^c*^ \dT/dC\ = oo. The 
singular behavior at v = Vs, C = C^, corresponds to the two peaks of Fig. 3a. 

The following is an explicit calculation of the value of ^ reached asymptotically for type 
II orbits (the calculation is similar for type I orbits). Expanding vq as Vq = — Svq = 
— I/xq — 6vq, we find from Eq. (fT8|) that the ffist order variations cancel exactly. The second 
order variations yield: 

v,{C)c^--il + \Ce'xl-l\'/') (B7) 

Xo 

(a similar calculation for type I events yields that the separatrix intersects the f — axis at the 
positive f— value f"^ = a^/xQ, = 0.278464... being the root of + Ina^ + 1 = 0; for type 
I curves the upper intersection with the v— axis is found through first variations of Eq. fllSj) . 
namely w7(C) = (I/xq) (0.278464 + 0. 108906 jCe^x^ _ i|)). 

The asymptotic behavior of the integral T(C) with respect to C can now be found by 
isolating the singularity of the integrand at f = fo(C) 

(-<g-2io^o(i _ 2xoA + 2xlA'^ + ...)dA 
{c)-vo(c)\ ^/i2xoCe-^^ovo + 2vo)A + {2Cxle-^^ovo _ i)A2 + ... 

where A = 1^; — Vq{C)\. The exact value of the lower limit used in this integral does not 
really matter in the calculation of dT/dC, since the leading contribution to T(C) comes from 
the parts of the orbits close to the X-point, i.e. for A small; the lower limit can actually 
be substituted by a value A 

max ~ 1/ Xq cusuring that the truncated expansion in the square 



2 ^ 
dC 



T{C) = 2 
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root is a sufficient approximation. On the other hand, it is necessary to retain O(A^) terms 
in the expansion within the square root of Eq. flB8|) . because the 0(A) term becomes very 
small as C tends to C^^ while the second order term is always of order unity. The upper 
limit of the integral fIBSp yields then logarithmic terms: 



T{C) ^ 2Ce-2^°'^« (Ji - 2x0/2 + 2xlh) (B9) 
with h = ln[8{xoCe-^^°'"> + Vq)] (2Cx^e-2*o^o _ iyi/2 

h = (3/2)(xoCe-2*«"o +t;o)'/n[8(xoCe-2*»"o +t;o)](2Ci;2e-2-o«o _ i)-5/2 _ 

As C ^ Cx, {xoCe~'^^°^° + Vq) — > 0, Ji becomes singular while I2 and Is are finite. Taking 
into account that (7e~^^°''° = Vq, we then find, to the leading order. 



or, using Eq.( lB6p and vq = —1 / xq — 5vq, 

i 1 



^0 

However, in view of Eqs. (]B7|) . (IBSP and flB4p we have Xo'^'^o °^ ^"^i; ^^^^ 

e 1 



^0 io^vi 

that is we obtain the power-law estimate of Eq. (l2T|) . 



+ ... (B12) 
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